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Multi-user diversity is considered when the number of users in the system is random. The complete 
monotonicity of the error rate as a function of the (deterministic) number of users is established and 
it is proved that randomization of the number of users always leads to deterioration of average system 
performance at any average SNR. Further, using stochastic ordering theory, a framework for comparison 
of system performance for different user distributions is provided. For Poisson distributed users, the 
difference in error rate of the random and deterministic number of users cases is shown to asymptotically 
^ approach zero as the average number of users goes to infinity for any fixed average SNR. In contrast, 

in 

for a finite average number of users and high SNR, it is found that randomization of the number of 
j— ( users deteriorates performance significantly, and the diversity order under fading is dominated by the 

smallest possible number of users. For Poisson distributed users communicating over Rayleigh faded 
channels, further closed-form results are provided for average error rate, and the asymptotic scaling 
law for ergodic capacity is also provided. Simulation results are provided to corroborate our analytical 
^ findings. 
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I. Introduction 

Point to point diversity combining schemes aim to mitigate the effects of fading in a wireless 
channel. In contrast, for multi-user systems another form of diversity termed multi-user diversity 
(MUD) is available, which thrives on the randomness of the user fading channels [1]. The key 
idea is to provide channel access to the user with the best channel at any instant of time. This 
has been shown to be optimal for both uplink [2] and downlink [3] scenarios. 

In the literature, MUD has been studied for the case of deterministic number of users only. 
Since the number of users is randomly varying in practice, it is of interest to consider MUD for 
this case as well. For example, cell phone users have longer voice calls while channel access 
for data communication is very short [4]. The probability of a cell phone user requesting data 
communication is very low, and bursty data requests such as stocks, weather and email lead 
to very short channel access times. This suggests that the number of users actively contending 
for channel access across time is random. Additionally, schemes in which a user is allowed 
to feedback its channel estimate to request channel access, when it is larger than a predefined 
threshold [5]-[7], also lead to a random number of users. Even in common scenarios where the 
fluctuations in the number of users is slower than the rapidity of channel fading, averaging error 
rates, or ergodic capacity, with respect to the user distribution results in meaningful system-level 
performance measures. 

In this paper, we analyze the performance of MUD systems with random number of users 
for the first time in the literature. In Section II the instantaneous SNR distribution of the best 
user chosen from a random set of users is derived for arbitrary fading and user distributions, 
and the mathematical preliminaries are presented. In Section III, the error rate averaged across 
fading with deterministic number of users is shown to be a completely monotonic function of 
the number of users N. Further, we also prove that the ergodic capacity of a MUD system with a 
deterministic number of users has a completely monotone derivative with respect to the number 
of users. These structural results of performance for a deterministic number of users are then used 
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to prove facts about the random number of users case. The first of these is that randomization 
of the number of users results in the deterioration of average performance measured in terms of 
either error rate or ergodic capacity, by using Jensen's inequality. In Section IV we introduce a 
framework in which different user distributions can be compared through the so called Laplace 
transform partial ordering of the number of users which is a particular stochastic order [8]. 
In Section V we derive the diversity order of a MUD system with random number of users 
and show that it is determined by the minimum possible number of users. In Section VI-A, 
expressions for outage for any fading distribution with Poisson user distribution are derived. For 
when the user distribution is Poisson distributed, Jensen's inequality for error rate is proved to 
be asymptotically tight in the average number of users in Section VI-B. For the special case 
when the number of users is Poisson distributed and when the user channel is Rayleigh faded, 
a closed-form expression for the error rate is derived in Section VI-C. The scaling of ergodic 
capacity with the average number of users is also provided for this case in Section VI-D. Poisson 
user distribution without allowing the number of users to be zero is considered in Section VI-E. 
Section VII corroborates our analytical results with simulations and Section VIII concludes the 
paper. 

Here are some remarks on notations used in this work. Asymptotic equivalence r(x) ~ g(x) 
as x — > a means that \im x ^ a r(x)/g(x) = 1, and r(x) = 0(g(x)) as x — > a means that 
limsup a ,_ 5>a \r(x)/g(x)\ < oo. In this paper, we consider a = or a = oo. Pr[-] is the probability 
of an event, E^[/(A')] is the expectation of the function /(•) over the distribution of the random 
variable X, and log(-) is logarithm to base e. 

II. System Model and Mathematical Preliminaries 

We consider an uplink MUD system with one base station (BS) and multiple users. Without 
loss of generality, both the BS and the users are assumed to have a single antenna. The received 
signal at the BS from the n th user can be expressed as, 

Vn = y/phnXn + W n , n = 0, 1, . . . , M, (1) 
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where the number of users Af is assumed to be a random variable with a discrete non-negative 
integer distribution. In Section VI-E we address the implications on the performance metrics 
of allowing the probability of Af = to be positive, leading to possibly no users and no 
transmission. When addressing the deterministic number of users case, we will set Af = N, 
where N is a realization of the random variable Af. A homogeneous MUD system is assumed 
where the average received power at the BS, p, is identical across all users. The symbol h n 
denotes the channel coefficient, x n the transmitted symbol, and w n the additive white Gaussian 
noise (AWGN) corresponding to the n th user. The channel is assumed to satisfy E[|/i„| 2 ] = 1 for 
all n and to be independent and identically distributed (i.i.d.) across all users. The transmitted 
symbols satisfy E[|x n | 2 ] = I. 

The channel gain of the n th user at the BS, prior to selection, can be expressed as 7„ = \h n \ 2 , 
and the selected user has a channel gain denoted by 7* = |/i*| 2 , where \h*\ 2 = max n {|/i n | 2 }. Note 
that 7* is a random variable that depends on the random variables Af, and \h n \ 2 , n = 0, 1, . . . , Af. 

Define F Jn (x) as the cumulative distribution function (CDF) of the channel gain of the n th 
user 7„. Since the fading channels across all users are assumed to be i.i.d., we drop the index 
n and define F 7 (x) := F ln (x). Recalling that the total number of users Af is a random variable, 
the CDF of the channel gain of the selected user, conditioned on Af = N, can be written as: 



where the N th power is obtained due to the i.i.d. assumption of the N user channels. The CDF 
of the channel gain of the best user selected from a random set of users can be obtained by 
averaging (2) with respect to the distribution of Af: 



where Uj^(t) = YlkLo ^ r [N = k]t h , < t < 1, is the probability generating function (PGF) of 
random variable Af. From (3) it can be seen that for any fading channel distribution and any 



F r (x\Af = N)=Fff(x), 



(2) 




(3) 
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non-negative integer distribution on the number of users, the CDF of the best user's channel 
gain at the BS can be easily obtained. 

We now survey some mathematical preliminaries that will be useful throughout. A function 
t(x) : IR + — > R is completely monotonic {cm.) if its derivatives alternate in sign [8], i.e., 

(-l) fc ^^>0, Vx, k = 0,1,2,..., (4) 

where d°T(x)/dx° = t(x) by definition. Due to a celebrated theorem by Bernstein [8], an 
equivalent definition for cm. is that it is a positive mixture of decaying exponentials. In other 
words, we have the Bernstein's representation t(x) = J °° e~ sx dip(s) for some nondecreasing 
function ip(s). In this paper, we are sometimes interested in cm. functions on integers, which 
are nothing but sequences obtained by sampling cm. functions as defined by (4). We are 
also interested in functions whose first-order derivatives satisfy (4), which are said to have 
a completely monotone derivative (c.m.d.). Even when the variable x is naturally an integer 
(such as the number of users), we will sometimes treat it as a real number, since we will be 
primarily interested in the asymptotic properties of t(x). 

A function t/j(s) is regularly varying with exponent /i ^ at s = oo if it can be expressed 
as ip(s) = s^l(s) where l(s) is slowly varying and by definition satisfies Hindoo l(ns)/l(s) = 1 
for k > 0. Regular (slow) variation of ip(s) at s = is equivalent to regular (slow) variation 
of ^(1/s) at oo. Intuitively, regular variation captures polynomial-like behavior near the origin 
or at infinity. The Tauberian theorem for Laplace transforms, whose proof can be found in [9], 
applies to cm. functions and states that t(x) is regularly varying at x = oo if and only if ip(s) 
is regularly varying at s = 0: 

Theorem 1: If a nondecreasing function ip(s) > defined on s > has Laplace transform 
t(x) = J °° e~ sx dil)(s) for x > 0, and l(s) is slowly varying at s = (or s = oo), the relations 
ip(s) ~ s^l(s) as s — > (or s — > oo) and t(x) ~ Y{ji + l)x~H{x~ l ) as x — > oo (or x — > 0) 
imply each other, where //Gl 
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One useful property given in [10, p. 27] is that 

/ s»l(s)ds ~ l(t) [ s^ds = -^—t^ +1 l(t) (5) 
Jo Jo + 1 

as t — > 0, with fj, > — 1 for l(s) slowly varying at s = 0. 

In this paper, we are interested in studying average error rates, and capacities averaged across 
both the channel distribution, and the number of users. The expression P e (p, N) represents the 
error rate of a MUD system with a deterministic number of users N, that is averaged with 
respect to the distribution of the fading channel. The expression E^- [P e (p,A/")] represents the 
average error rate of a MUD system with a random number of users, which is averaged with 
respect to the distribution of the number of users and the fading channels. 

III. Properties of the Average Error Rate and Ergodic Capacity 
A. Average Error Rate 

In this section, we first prove that the average error rate of a MUD system, with a deterministic 
number of users N, is a cm. function of N, under general conditions. This will be used to infer 
about the behavior of the average error rate when a random number of users is considered, in 
Section IV. 

The error rate of a MUD system with a deterministic number of users N and average SNR 
p is given by, 

roc 

P e (p,iV)= / P e (px)dF^(x) (6) 
Jo 

where P e (p:r) is the instantaneous error rate over an AWGN channel for an instantaneous SNR px 
of the best user. Often, the instantaneous error rate is assumed to have the form P e (pa;) = ae~ vpx 
or P e (px) = aQ(y/rfpx), where a and rj can be chosen to capture different modulations [11]. To 
represent (6) in terms of the CDF F y (x), rather than the probability density function (PDF), we 
left it as a Stieltjes integral [12] even though it can also be expressed in terms of the PDF / 7 (x) 
using dF^(x) = NF^~ 1 (x)f 1 (x). In what follows, we will study the sequence P e (p, N) as a 
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function of the integer variable N. Since we are ultimately interested in the asymptotic behavior 
of this sequence, we will also consider (6) with N being a real number. 

We begin by proving that P e (p, N) is a cm. function of iV not just for P e (px) in the forms 
of exponential function and Q function, but for any instantaneous error rate function. In other 
words, we only assume ¥ e (px) is decreasing in x for any p > 0. Defining B(x) = —dP e (x)/dx, 
after integrating (6) by parts, the k th derivative of P e (p, N) can be written as, 



Since P e (pa;) is decreasing, and log (F 1 (x)) < we see that (7) satisfies the definition in (4). In 
particular, P e (p, N) being a cm. function means that (7) is negative for k = 1 and positive for 
k = 2, and consequently P e (p, N) is a convex decreasing function of N. For when the number 
of users in the system is random, by applying Jensen's inequality for convex functions, we have, 



where A := E[Af]. Therefore, randomization of the number of users always deteriorates the 
average error rate performance of a MUD system. 

To establish the complete monotonicity of P e (p, N) as a function of N, we only used the fact 
that the instantaneous error rate P e (p^) in (6) is a decreasing function of x for p > 0, which 
always holds. This cm. property will be used to stochastically order user distributions in Section 
IV. 

B. Ergodic Capacity 

The ergodic capacity for the deterministic number of users system can be expressed as, 



where we use integration by parts, and assume that F 1 (x) satisfies lim^oo log(l + px)(l — 
F™{x)) = 0, for all N > 0. 




(7) 



E M [P e (p,Af)] >P e (p,A), 



(8) 




(9) 
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It can be seen that C(p, N) has a completely monotonic derivative since, 



C(p,N) 




(x) [log(F,(x))} 



k+l 



dN k + l 



dx. 



(10) 



1 + px 



alternates in sign as k is incremented. This establishes that C(p, N) has a completely monotonic 
derivative, provided that the fading distribution satisfies the mild assumption lim^o log(l + 
px)(l — F^(x)) = for all N > 0, as assumed after (9). This assumption holds for all 
distributions with exponential or power law tails, which is the case for all fading distributions in 
wireless communications. Using (10) with k — 0, 1, it is seen that C(p, N) is concave increasing 
function of N. Applying Jensen's inequality for concave functions, we have 



Therefore, similar to the error rate metric, randomization of N will always hurt the average 
ergodic capacity of a MUD system. The c.m.d. property of C(p, N) will be used in the following 
section discussing the stochastic Laplace transform ordering of user distributions. 



In this section we introduce Laplace transform (LT) ordering, a tool to compare the effect 
that different user distributions has on the error rate, and ergodic capacity averaged across user 
and channel distributions. Stochastic ordering of random variables, of which LT ordering is a 
special case, is a branch of probability theory and statistics which deals with binary relations 
between random variables [8], [13]. 

Let X and y be non-negative random variables. X is said to be less than y in the LT order 
(written X < Lt y), if E [e~ s *] > E [e~ sy ] for all s > 0. An important theorem found in [8], 
and [13] is given next: 



Theorem 2: Let X and y be two random variables. If X < Lt y, then, E [^(A")] > E [^(3^)] 
for all cm. functions ip(-), provided the expectation exists. Moreover, the reverse inequality 



E M [C(p,Af)} <C(p,X). 



(11) 



IV. Laplace Transform Ordering of User Distributions 
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E['0(Af)] < E[^(y)] holds for all ip(-) with a completely monotone derivative, provided the 
expectation exists. 

It follows that if two user distributions satisfy Mi <u Mi, then for all average SNR p 

E M [Pe(p,M)] > E^ 2 [? e {p,M 2 )] , 

E M p(p,M)] < E^ p(p,JV 2 )] . (12) 

To rephrase (12), if the number of users M is from a distribution that is dominated by another 
distribution in the Laplace transform sense, then both the average error rate and capacity are 
respectively ordered at all average SNR p. 

The LT ordering of discrete random variables can also be expressed in terms of the ordering 
of their PGFs. By defining t := e~ s , one can rewrite E [e _s *] > E [e^ 5 ^] for s > as 
E [t x ] > E [t y ] for < t < 1, which is the same as U x {t) > U y (t), < t < 1, where we 
recall that U x (t) = E[t x ] represents the PGF of the discrete random variable X. 

To provide examples of random variables that are LT ordered, consider Poisson random 
variables X and y with means A and p respectively, such that A < p. It is straightforward to show 
that for this case e x ^ = E [**] > E [t y ] = e^^K for < t < 1, implying that X < Lt y. If 
X and y are geometric distributed with probability of success on each trial pi and p 2 respectively, 
such that pi < p 2 , then y < Lt X since p 2 /(l - (1 -p 2 )t) = E [t y ] < E [t x ] = pi/(l - (1 -pi)t) 
for < t < 1. Similarly, for X being Poisson distributed with parameter A and y being geometric 
distributed with parameter p = 1/(1 + A), so that E[X] = E[y] = A, it can be once again shown 
that E [t y ] < E [t x ] for < t < 1, establishing y < Lt X. From this latter result, one can 
conclude that Poisson offers a better user distribution than geometric distribution for a fixed 
average number of users at all average SNR p, from both error rate and capacity points of view. 

V. High SNR Analysis and Diversity Order 

In this section we analyze the average error rate at high SNR under general assumptions on 
the user distribution and fading channel distribution. In the following, we will assume that the 
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fading distribution F 1 (x) is regularly varying with exponent d at x = (typically F 7 (x) = 0(x d ) 
as x — > 0), which is true for many fading distributions including Rayleigh (d = 1), Nakagami-m 
(d = m) and Ricean (d = 1) [14]. 

Theorem 3: Let F 7 (x) be regularly varying at x = with exponent d > 0, and Af G {/c , k + 
1, ...} represent the range of the number of users random variable. Then the high-SNR asymptotic 
average error rate is given by 

[P e (p, Af)] ~ Pr [A/" = fc„] CiF^^p" 1 ) (13) 

as p — >■ oo, where Ci and C 2 are constants given by C\ = aT(k d + 1), C 2 = when 
P e (px) = ae-^ px and d = aT(k d + l/2)/(2y/n), C 2 = 2^ when P e (px) = aQ(^/rjpx). 

Proof: See Appendix A. ■ 
From (13) it is straightforward to show that the diversity order of the MUD system with 
random number of users is given by k d. This follows from (13) and the regular variation 
assumption on F 7 (x), so that F*°(C 2 p~ 1 ) = O(p~ kod ) as p — > oo. 

VI. Poisson Distributed M 

Consider a MUD system which contains a large number of users. Suppose each user is 
active with a small probability. In such a system, as the number of users increases, the user 
distribution will approximate the Poisson. In this section we analyze the system when Af is 
Poisson distributed with parameter A. 

A. Outage Probability and Its Asymptotic Behavior for Large A 

When Af is Poisson distributed with parameter A, the probability of outage with a threshold 
x can be expressed as, 

F r (x) = Pr[ 7 * <x] = Ve- A -Fji) = e -W-M*)) i[ x > ] (14) 

fc =o k - 
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where /[•] is the indicator function. Equation (14) implies that as the average number of users 
increases, the outage probability decreases for any distribution on 7„, the channel gain of the 
user fading channel. 

In what follows, we show that for large A the outage behavior is dependent on F 7 (x) only 
through its tail behavior. In fact, it is possible to show that there exist normalizing and shift 
functions a(A) and 6(A) such that the probability Pr [(7* — 6(A))/a(A) < x] for large A is 

lim F r (a(\)x + 6(A)) = exp(-e~ x ), -00 < x < 00 (15) 

which is known as the Gumbel distribution [15]. Using (14), sufficient and necessary conditions 
for (15) are clearly that lim^^ A(l — F 1 (a(X)x + b(X))) = e~ x . In [15, p. 300] it is shown 
that a(A) = [A/ 7 (6(A))] _1 , 6(A) = — 1/A), satisfy this condition for many distributions 

including Rayleigh, Nakagami-m and Ricean. Also, it can be seen that the asymptotic CDF of 
the Poisson number of users case for large A has the same form (Gumbel distribution) as the 
asymptotic CDF of the deterministic number of users case for large N. 

B. Average Error Rate 

In our outage analysis for the Poisson number of users case, we showed that the outage 
probability in (14) for large A approaches the Gumbel distribution, which is also the asymptotic 
distribution obtained for a deterministic number of users case. Therefore, even though we saw 
that randomization always deteriorates performance, for large average number of users it should 
approximately yield the same performance as the deterministic case. This amounts to the tightness 
of Jensen's inequality for the Poisson users case. 

We now provide sufficient conditions for Jensen's inequality involving P e (p, N) in (12) to 
be asymptotically tight in A. Recall that P e (p, N) is the error rate averaged over the channel 
distribution for deterministic number of users N. To this end, we use the results in [16, Theorem 
2.2] which were derived in a networking context for arbitrary cm. functions. 

Theorem 4: Let P e (p, iV) be cm. and regularly varying at N = 00 and consider the error 
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rate averaged across the channel and the users [P e (p, AO] , where J\f is a Poisson distributed 



random variable with mean A. Then, 



Eat [P e (p, AO] = P e (p, A) + O (P e (p, A)/A) 



(16) 



as A — > oo. 



Equation (16) shows that as A — > oo, the difference between the error rate averaged across 
the user distribution and the error rate evaluated at the average number of users vanishes as A 
tends to oo. This implies that for sufficiently large A the performance of the MUD systems with 
random number of users will be almost equal to the performance of the MUD systems with a 
deterministic number of users with the number of users equal to A. 

To apply Theorem 4 we require P e (p, N) to be cm. and regularly varying. We have already 
shown that P e (p, N) is always completely monotonic in N. Next, we provide the conditions 
under which P e (p, N) is a regularly varying function of N. Consider 



Jo 

where B(-) is defined as B(x) = —dP e (x)/dx. Now, setting u := — log(F 7 (x)), and integrating 
by substitution we have, 



where F~ 1 (x) is the inverse CDF and / 7 (x) is the PDF of 7„. We now establish the sufficient 
conditions for P e (p, N) to be a regularly varying function of N: 

Theorem 5: If P e (p, N) is cm. in N, a sufficient condition for it to be regularly varying at 
iV = oo is that, t{u) := p{B{pF- 1 {e- u ))e- u )/{f 1 (F^e - "))) is regularly varying at u = 0. 

Proof: By comparing the representation of P e (p, N) in (18) with the Bernstein's represen- 
tation of cm. functions discussed after (4), it can be seen that (18) can be represented as the 
Laplace transform of t(u). Using Theorem 1, the proof follows. ■ 




(17) 




(18) 
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Theorem 5 shows that for the conclusions of Theorem 4 to hold (i.e., Jensen's inequality to 
be asymptotically tight), the CDF of the single-user channel F 7 (x), and the error rate expression 
P e (px) have to jointly satisfy the regular variation condition given in Theorem 5. Next, we 
examine whether this condition holds for commonly assumed instantaneous error rates P e (px) 
with 7 n being exponentially distributed. For the case of ¥ e (px) = ae~ vpx , we have t(u) = 
ap(l — e~ u ) 1JP ~ 1 e~ u , which satisfies lim u ^. t(Ku)/t(u) = K np ~ 1 , therefore proving the regular 
variation of t(u) at 0. By using Theorem 1 this in turn proves regular variation of P e (p, N) 
at N = oo. Therefore P e (p, N) is both a cm. and a regularly varying function of N for this 
case. Consequently, when P e (px) = ae~ vpx and the fading is Rayleigh (i.e. channel gain is 
exponential), the difference in error rate performance of a MUD system with a random number 
of users averaged over the number of users distribution and of a deterministic number users 
approaches zero for sufficiently large A, as in Theorem 4. 

Consider now P e (px) = aQ(^frfpx), with 7„ being exponentially distributed. The error rate 
can be expressed as, 



/(i 2^2n Jo ^fx 



poo 

>Xp,N) = a / Q(^)dF?(x) 
Jo 



P e (p, N) = a Q {yfqpx) dF" (x) = —^== / = dx, (19) 



where the second equality is obtained by integration by parts. Once again, by setting u 
— log(l — e~ x ) we can rewrite (19) as, 



awrjp 



roo -u 

_ / exp (-Nu) (1 - e~ u ) r ' p,2 ~ l -. =du. (20) 

2V^Jo V- log(l-e-«) 

Thus we have t(u) = a^/rfp(l — e ~ u ) , ' /5 / 2_1 e _ "/ (2a/— 2tt log(l — e^ u )) and it can be shown that 
\im u ^,o t(K,u)/t(u) = k^p/ 2 ^ 1 , therefore once again proving that P e (p,iV) is both a cm. and a 
regularly varying function of N. Having verified the conditions of Theorem 5 for P e (px) = 
aQ(y/rfpx) with 7„ being exponentially distributed, we conclude the tightness of Jensen's 
inequality as suggested by Theorem 4. 
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C. A Special Case: Poisson distributed Af and Rayleigh Faded Channels 

In this section, we consider the case when the number of users Af is Poisson distributed and 
the user channels are Rayleigh faded. This practically relevant case will lead to closed form 
expressions. 

1 ) Distribution of Channel Gain: For this case the CDF of the channel gain of the best user 
using (14) is given by, 

F r (x) = exp (-Xe~ x ) I[x > 0]. (21) 

The channel gain of the best user in (21) is identical to a truncated Gumbel distribution, which 
was seen in its untruncated form in (15). Notice that for x = (21) yields e~ A > so F^*(x) 
has a jump at x = 0. The distribution in (21) is therefore of mixed type with a mass of e~ A at 
the origin and the rest of the distribution has the form of a truncated Gumbel distribution. 

2) Average Error Rate: Assuming the error rate has the form, P e (px) = ae~ vpx as mentioned 
in Section III-A, the average error rate can be expressed as, 

POO POO 

E M [P e (p,A0] = Aa / e-i^e-'e-^'dx + a S(x)e~ > 'e"^ 'dx . (22) 
Jo Jo 

Setting y = \e~ x and integrating by substitution, (22) can be expressed as, 

/•X 

[P e (p, A/")] = a J ' e~ y dy + ae~ x = aX'^rjp + 1, A) + ae~\ (23) 

where 7(5, x) is the lower incomplete gamma function [12]. It can be easily shown that a\~ r}P ^(r]p+ 
1, A) + ae~ x ~ a\~ np Y(r]p + 1) as A — > 00, indicating a power-law decay in the error rate as a 
function of the average number of users. 

D. Asymptotic Scaling of Capacity with A 

Next, we derive the asymptotic average capacity and the corresponding scaling laws with 
respect to A. 
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Theorem 6: For Poisson distributed Af with mean A and Rayleigh faded channels, as A — > oo, 
we have 

[C(p,Af)} = log(l +plog(A)) + 0(1/ V / MA)). (24) 
Proof: See Appendix B. ■ 
For a MUD system with deterministic number of users N, it has been shown in [1] that the 
ergodic capacity grows as loglog(iV). From Theorem 6 it is seen that for a MUD system with 
random number users, whose mean is A, the ergodic capacity grows as loglog(A). This implies 
that when average number of users A is equal to N of the deterministic number of users case, 
the ergodic capacity for both cases grow at the same rate. 

E. Zero Truncated Poisson User Distribution 

The CDF expression in (21) includes the case when N = 0, i.e., there are no users in the 
system. When there are no users, no data will be transmitted. In view of this, it is reasonable to 
drop the N = case and model the user distribution with the zero-truncated Poisson distribution 
which is given by Pr[A/* = k\Af > 0] = Pr[jV = k]/(l — Pr[jV = 0]), for any positive integer 
k. For zero-truncated Poisson distributed Af, N e {1,2, . . .} and mean A = A/(l — e~ A ) where 
A is the mean of the underlying Poisson random variable. The CDF of the channel gain of the 
best user can be expressed as, 

1 ~ P -A P W--e~ x ) _ 1 

FAx) ^E^fAV = e A_! • ( 2 5) 

k=l 

For when the average number of users A — > oo, it can be seen that A ~ A. Further, 

A(l— e~ x ) i 

lim FAx) = lim — — = e~ x{e ~ x \ (26) 

which is identical to the CDF for Poisson distributed user case in (21). This implies that for a 
large average number of users in the system, the outage, average error rate, and ergodic capacity 
performance of zero-truncated Poisson distributed user case will be identical to that of the Poisson 
distributed user case. 
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VII. Simulations 

An uplink MUD system where both the BS and users having a single antenna is considered. In 
this section, using Monte-Carlo simulations, the error rate, ergodic capacity and outage capacity 
are simulated to corroborate our analytical results. For all simulations considered, the Rayleigh 
fading is assumed. 

In Figure 1, assuming n/A QPSK modulation, the average bit error rate with deterministic TV 
is compared with the performance averaged across various user distributions. It is seen that the 
deterministic number of users system performs better than all the cases involving random number 
of users. The performance of the Poisson distributed users case comes close to the deterministic 
case as A increases, as predicted by Theorem 4. 

In Figure 2, the ergodic capacity is plotted against A for the random cases and N = A for 
the deterministic case. It is seen that the capacity of the deterministic number of users system 
is the highest while for all distributions of J\f, the capacity is worse, corroborating our result in 
Section III-B. 

In Section IV, we showed that Poisson distributed random variables and geometric distributed 
random variables are LT ordered, which also orders their respective average error rate and ergodic 
capacities when averaged across the respective user distributions. In Figures 3 and 4 it can be 
seen that both error rate and capacity follow their corresponding ordering at all average SNR p. 

In Figure 5 the bit error rate versus average SNR for different average number of users is 
shown. It can be seen that the analytical approximation of average error rate with Poisson number 
of users derived in (23) is within 1 dB of the Monte-Carlo simulation result. Following the result 
in (23) it can also be seen that the larger the value of A, the lower the error rate. 

In previous simulations we saw that increasing A leads to an improvement in performance, 
and for a fixed average SNR the performance of the system with Poisson distributed number 
of users approaches the performance of the system with deterministic number of users. Figure 
6 considers a zero-truncated Poisson distribution and illustrates that at high average SNR, the 
diversity order is k d = 1, verifying Theorem 3. This leads us to conclude that for low SNR's 
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but sufficiently large A the performance of the random number of users is nearly identical to 
that of the deterministic case. However for high SNR's, the performance of the random number 
of users case is significantly worse due to the loss in diversity order. 

VIII. Conclusions 

Multi-user diversity (MUD) is analyzed for when the number of users in the system is random. 
The error rate of MUD systems is proved to be a completely monotone function of the number 
of users in the system, which also implies convexity. Further, ergodic capacity is shown to have a 
completely monotone derivative with respect to the number of users. Using Jensen's inequality, 
it is shown that the average error rate and ergodic capacity averaged across fading and the 
number of users will always perform inferior to the corresponding performance of a system with 
deterministic number of users. Further, we provide a method to compare the performance of the 
system for different user distributions, using a specific stochastic ordering based on the Laplace 
transform of user distributions. 

Importantly, for the MUD system with random number of users, it is shown that the diversity 
order is defined by the minimum of the range of realizations of the number of users. When 
the number of users are Poisson distributed, for any user channel fading distribution, outage 
probability is shown to converge to the truncated Gumbel CDF, similar to the case of the 
deterministic number of users system. Further, it is proved that the difference between the error 
rate performance of the Poisson number of users system and the deterministic number of users 
case goes to zero like O (P e (p, A)/A) asymptotically in the average number of users. As a special 
case, when the user fading channels are Rayleigh distributed, a closed-form error rate expression 
is provided. Also, the asymptotic scaling law of ergodic capacity is analyzed and shown to be 
approximately log (1 + plog(A)) for large A. Finally, zero-truncated Poisson number of users 
case is shown to not affect our main conclusions for the common scenario where the number 
of users is always positive. 
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Appendix A. Proof of Theorem 3 

Since F 7 (x) is regularly varying, it must be in the form F 7 (x) = x d l(x) where l(x) is slowly 
varying at 0. For a system with k users, the CDF of the channel gain becomes F k (x) = x kd l k (x). 
It is easy to verify that l k {x) is slowly varying at 0. Therefore, given j Q l dF k {x) ~ t kd l k (t) and 
P e (px) = ae~ vpx , it follows based on Theorem 1 that 

POD 

P e (p, k) — a e~^ x dF k {x) ~ «r(A;rf + l)F fc (jrV _1 ) (27) 

as p — y oo. For the case P e (pa;) = aQ(^frfpx), the asymptotic expression of P e (p,/c) can be 
derived similarly as follows. Using integration by parts we obtain 

P e (p, k) — a Q(^)dF k (x) = -^f / -^F k (x)dx (28) 
Jo ly Itx Jo v x 

Based on (5) we have 

jT -^F 7 fe (x)dx ~ / fc (t) jT x ^ 2 dx = kd l 1/2 tkd+1/2lk ^ (29) 

as t — y 0. Then the asymptotic average error rate given by (28) becomes 



(30) 



aT(kd+ 1/2) (r] P \-kd aT(kd+ 1/2) . _ x _ lN 



2^ V 2 / v ' r y 2^ 
as p — >■ oo, where the asymptotic equality is based on Theorem 1 and the second equality is 
based on F 7 (x) = x d l(x). It can be seen that both (27) and (30) have the form P e (p, k) ~ 
C 1 F k (C 2 p~ 1 ) with C\, C 2 being constants. Consequently, as p — y oo, the dominant term in the 
average error rate E^y [P e (p, A/")] = YlkLk ^ r \^ = ^1 Pe(p 5 is the term with k = k , which 
deceases slower than any other term. Using dominated convergence theorem, we can easily 
determine lim^^ E^ [P e (p, A/")] /P e (p, k ) = Pr [Af = k ] by exchanging limit and summation. 
Consequently, the ratio between the left hand side and the right hand side of (13) goes to 1 as 
p — y oo, and we thus have the asymptotic average error rate given by (13), with C\, C 2 given 
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as in the theorem. 

Appendix B. Proof of Theorem 6 

For Poisson distributed Af and Rayleigh faded channels, using integration by parts and the 
CDF in (21), the capacity of the system can be written as, 

E X [C(p,jV)] = P / -r- dx. (31) 

Jo i + px 

Defining y := e~ x and integrating by substitution, 



= (^ )/X 1 - C ' A " d y + rf (32) 

7o 1 -Plog(y) V2// A/ma)/a1/(1 -Plog(y)) 



For the first term after the second equality in (32), we have 

ry^W/x 1 _ e -xy / p \ rv^A Xy /p\ 

< Jo I-Plog(y) \y) y< J l+plog(A)-( P /2)log(log(A)) \y) 

= p V / log(A) 

l + plog(A)-(p/2)log(log(A))' 

(33) 

by replacing the numerator of the integrand with its upper bound and the denominator of the 
integrand with its lower limit. It can be seen that the upper bound after the equality in (33) 
yields 0(l/-\/log(A)) and has limit as A — > oo, implying that the first term should have limit 
0. The second term in (32) has the bounds given by, 



A/waI/a y(l-plog(y)) A a 



'Vbgw/A y(i-pi°g(y)) i^/bgw/A^ 1 -piog(y)) ivw/A -P lo s(y)) 

(34) 

in which the lower and upper bounds are obtained by bounding the numerator, and they turn 

out to be (l - e ~ v 71 ^) j l g(l + p log(A) - (p/2) log(log(A))) and (1 - e~ x ) log(l + p log(A) - 
(p/2) log(log(A))) respectively. Also, it can be verified that log(l+p log(A) — (p/2) log(log(A))) = 
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log(l+plog(A)) + 0(log(log(A))/log(A)) as A ->■ oo, and lim^ e'V 71 ^) log(l + plog(A)) = 
liniA^oo e _A log(l + plog(A)) = 0. Therefore, for a fixed p, and as A — > oo we can express (32) 
as (24) considering the fact that log(log(A))/log(A) decays faster than l/y/\og(\), completing 
the proof. 
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Fig. 1. Error rate vs. A: Rayleigh Fading Channel, average SNR = 6 dB 
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Fig. 2. Capacity vs. A: Rayleigh Fading Channel, average SNR = lOdB 
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Fig. 3. Error rate vs. average SNR: Rayleigh Fading Channel 




Fig. 4. Capacity vs. average SNR: Rayleigh Fading Channel 
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Rayleigh faded channel with parameter 1 
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Fig. 5. Error rate vs. average SNR: Poisson Users and Rayleigh Fading Channel 
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Fig. 6. Diversity Analysis: Poisson Users and Rayleigh Fading Channel 
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